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Abstract 

The  dielectronic  recombination  rate  coefficients  have  been 
calculated  for  the  various  ionization  stages  of  silicon.  Account 
has  been  taken  of  all  stabilizing  radiative  transitions  and  all 
autoionization  processes  which  involve  a single-electron  electric- 
dipole  transition  of  the  recombining  ion  core.  For  certain  ions 
the  dielectronic  recombination  rates,  although  still  larger  than 
the  direct  radiative  recombination  rates,  are  found  to  be  substanti- 
ally reduced  when  account  is  taken  of  the  effects  of  a previously 
neglected  autoionization  process  in  which  the  excited  recombining  ion 
core  undergoes  a &n  = o transition  to  a lower  excited  state.  The 
temperatures  at  which  these  ions  have  their  maximum  abundance  in 
corona  equilibrium  are  significantly  reduced  when  use  is  made  of 
the  new  dielectronic  recombination  rates.  Calculations  are  also 
presented  for  the  total  rates  of  radiative  energy  loss  from 
isothermal  steady-state  plasmas  due  to  the  line  and  continuum 
emission  of  silicon  ions. 
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I.  Introduction 

An  important  method  of  investigating  the  physical  properties  of 
high-temperature  laboratory  and  astrophysical  plasmas  is  the  analysis 
of  optically-thin  spectral  lines  emitted  by  multiply-charged  ions.  In 
addition,  the  energy  loss  due  to  the  emitted  far  ultraviolet  and  x-ray 
radiation  may  have  a significant  influence  on  the  dynamical  behavior  of 
the  plasma. 

At  low  densities,  where  practically  all  excited  ions  cascade  to 
their  ground  states  in  times  that  are  short  compared  with  the  total 
electron-ion  collision  time,  the  spectral  line  intensities  can  be  pre- 
dicted on  the  basis  of  the  corona  model  (McWhirter  I965) . The  steady- 
state  abundance  of  the  various  ionization  stages  of  a given  element  can 
be  determined  as  functions  of  the  local  electron  temperature  by  assuming 
that  electron  impact  ionization  from  the  ground  states  is  balanced  by 
direct  radiative  and  dielectronic  recombination.  Photoionization  and 
stimulated  radiative  recombination  processes  are  neglected  because  the 


plasma  is  assumed  to  be  transparent  to  the  emitted  radiation.  The  spectral 


line  intensities  can  then  be  calculated  from  the  corona  ionization 


equilibrium  distributions  and  the  electron  impact  excitation  rates. 

Dielectronic  recombination  has  been  shown  (Burgess  1900  to  be  the 
dominant  recombination  process  for  nonhydrogenic  ions  in  low-density 
high-temperature  plasmas.  In  addition,  Gabriel  and  Jordan  (I972)  have 
emphasized  the  importance  of  determining  the  contribution  to  the  resonance 
line  intensity  due  to  spectroscopically  unresolvable  dielectronic  recom- 
bination satellites.  Burgess  (I965)  has  obtained  a simple  expression 
for  estimating  the  dielectronic  recombination  rate  coefficients  for 


■ ~ — 


t—  *■  3 — ' ^ r. 3^ — *1  l.f 


A 


-4- 

low-density  plasmas.  The  Burgess  formula  has  been  employed  in  the 
calculation  of  the  corona  ionization  equilibrium  distriubtions  (Jordan 
I969)  and  in  the  evaluation  of  the  radiative  energy  loss  rates  due  to 
collisionally-induced  line  emission  (Cox  and  Tucker  1969),  and  the 
importance  of  dielectronic  recombination  has  been  firmly  established. 

In  this  investigation  a more  general  theory  (Jacobs,  Davis,  Kepple, 
and  Blaha  1977),  referred  to  subsequently  as  JDKB,  has  been  applied  in 
systematic  and  detailed  calculations  for  the  dielectronic  recombination 
rates  of  the  various  ions  of  silicon.  The  Burgess  formula  is  found  to 
overestimate  the  dielectronic  recombination  rates  for  certain  ionization 
stages  because  of  the  neglect  of  the  effects  of  an  autoionization  process, 
not  previously  discussed,  in  which  the  excited  recombining  ion  core  under- 
goes a in  = o single-electron  electric-dipole  transition  to  a lower  excited 
state.  The  new  dielectronic  recombination  rates  are  used  to  calculate 
the  corona  ionization  equilibrium  distributions.  Finally,  these  distri- 
butions are  employed,  in  conjunction  with  reliable  distorted  wave  results 
for  the  electron  impact  excitation  rates  (Davis  et  al.  1977) > in  a new 
evaluation  of  the  radiative  energy  loss  rates. 


II.  Theory  of  Dielectronic  Recombination 

+Z 

Dielectronic  recombination  of  the  ion  X with  residual  charge  Z is 


the  result  of  a radiationless  capture  into  a doubly-excited  state  j,  n£ 
_ V+(Z-1) 
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followed  by  a stabilizing  radiative  transition  to  a state  k,  n£  which 


lies  below  the  ionization  threshold 


x+(z-i)  (j>  ni)  ^ x+(z-i)  (k>  ^ + h ^ 


In  the  corona  model,  the  recombining  ions  are  assumed  to  be  initially 
in  their  ground  states  (i  = g) . In  addition,  the  most  important  stabi- 
lizing radiative  transitions  will  return  the  recombining  ion  core  to  its 
ground  state  configuration  (k  = g) . In  this  investigation  the  many- 
electron  states  i,  j,  and  k will  be  specified  by  the  effective  principal 
and  angular  momentum  quantum  numbers  of  the  active  electron. 

For  a Maxwellian  electron  velocity  distribution  the  dielectronic 
recombination  rate  coefficient  into  the  final  state  k,  nl  can  be  ex- 
pressed in  the  form  (Shore  I969) 


— o3  _ l"2-  rc  l\r  T ^2 


ad  (Z,  i Z-lj  k,  m ) = 23  a^  v * (E^  Tg) 


x v-*  g(Z-l;  j , vJS)  ^a( j , n-L  -♦  i,E A^( j ,n-L  -♦  k,nt,) 
, 2 g(Z, i)  A(Z-1;  j,ot) 


X exp 


E(Z,i)  - I 
kBTe 


i)  - E(Z-1 ; j,nt) 


The  statistical  weights  associated  with  the  energy  levels  E (Z,i)  and 

E ( Z— 1 ; j,ivl)  are  denoted  by  g (Z,i)  and  g (Z-l;  j,n£)  and  E„  = 13.6  ev. 

For  Z S 20  dielectronic  recombination  occurs  primarily  through 

transitions  involving  highly  excited  nl  states  of  the  outer  electron, 

especially  when  the  recombining  ion  core  undergoes  a An  =0  transition. 

Consequently,  the  autoionization  rate  A (j,n£  -+  i,  6 ),  in  terms  of  which 

3 1 


m 


the  radiationless  capture  rate  has  been  expressed,  can  be  obtained  from 
the  threshold  values  of  the  partial  wave  electron  impact  excitation 
cross  sections  a (i,^  t.  -*  j,  t^)  by  means  of  the  quantum  defect  theory 
relationship  (Seaton  I969) 

A ( j , nt  -4  i,e  ) = fH  __gL2 

h Tr2  n3  a2 
o 


[ E(2»j>  ~ E(Z,i)  1 g(Z.i) 

L Eh  2(21  + 1)  g(Z,j) 


22  ■*  j>  l) 


j = O 


In  addition,  the  spontaneous  radiative  transition  rate  A^Cj,  nt  -♦  k,  nt) 
can  be  calculated  by  neglecting  the  influence  of  the  outer  nt  electron. 

The  essential  innovation  in  the  present  investigation  is  the  correct 
determination  of  the  total  decay  rate  A(Z-1;  j,  nt) . A careful  exami- 
nation of  the  general  theory  has  revealed  (JDKB)  that  the  total  decay 
rate  of  the  doubly-excited  state  j,  nt  must  be  determined  by  taking  into 
consideration  the  autoionization  rates  A (j,  nt  -»  k,e,  ) into  all  lower 
states  k of  the  recombining  ion  and  the  spontaneous  radiative  transition 


rates  A (j,  nt  -»  k,  nt)  and  a (j,  nt  -»  k,  8 , ) into  all  lower 
r r k 

states  k,  nt  and  continuum  states  k,  e,  . 

k 


bound 


A(Z-1;  j,  n-1)  = 


£ ka 

‘ < j L 


Ar(j,  nl  -*  k,  nl)  + j dek  A^j,  n-t  -»  k,  ek)  . 


In  previous  calculations  it  has  been  assumed  that  the  doubly-excited 
states  which  give  the  most  important  contributions  to  the  total  dielec- 
tronic  recombination  rate  undergo  autoionization  with  the  greatest 
probability  into  the  ground  state  of  the  recombining  ion.  While  the 
most  probable  radiative  decay  process  usually  terminates  in  the  ground 
state,  the  corresponding  assumption  for  autoionization  is  invalid  for  the 
most  important  doubly-excited  states  of  certain  ions.  In  the  boron-like 
through  the  neon-like  ions,  for  example,  the  Jd,  nl  doubly-excited  states 
are  found  to  autoionize  preferentially  into  the  lower  excited  3p  state 
rather  than  into  the  2p  ground  state,  from  which  radiationless  capture 
is  assumed  to  occur.  This  is  illustrated  in  Figure  1 by  using  long  and 
short  arrows  to  indicate  the  alternative  autoionization  channels.  Because 
capture  from  the  3p  excited  state  is  negligible  at  low  densities,  as  in- 
dicated by  the  absence  of  the  corresponding  arrow  in  Figure  1,  the  in- 
clusion of  the  additional  3d  -»  3P  autoionization  rate  must  lead  to  a 
reduction  of  the  Jd  -*  2p  contribution  to  the  dielectronic  recombination 
rate. 

Previous  investigations  of  dielectronic  recombination  contain  no 

o 

discussion  of  the  role  of  radiative  Auger  transitions  (Aberg  1971)  in 
the  determination  of  the  total  decay  widths  A(Z-1;  j,  nt) . Using 
quantum  defect  theory  (Seaton  I969)  the  total  contribution  from  radiative 
Auger  transitions  of  the  type 


-8- 


X+(Z-1)  ^ ^ X+(Z)  (k>  + e-  (e^k)  + * m (6) 

is  found  to  satisfy  the  inequality 

| d*k  Ar(j,  nt  -♦  k,  ek)  < 2|f_  Ar(j-*k).  (7) 

It  follows  that  for  large  n-values  the  radiative  Auger  contribution  to 
the  total  decay  width  can  be  neglected  provided  Z is  not  too  large. 

The  basic  approximation  in  the  theory  of  dielectronic  recombination 
is  that  the  total  spontaneous  radiative  capture  rate  coefficient 
a (Z,  i -»  Z-l;  k,  n-t)  into  the  final  state  k,  n-C,  can  be  obtained  simply 
by  adding  to  the  dielectronic  contribution  given  by  equation  (3)  the 
direct  radiative  recombination  rate  coefficient  describing  the  free- 
bound  transition  of  a single-electron.  A formula  for  estimating  the 
direct  radiative  recombination  rate,  which  takes  into  account  the  presence 
of  equivalent  electrons  in  the  outer  subshell  of  the  recombining  ion,  has 
been  given  previously  (JDKB).  The  interference  between  the  direct 
radiative  and  dielectronic  recombination  amplitudes,  which  occurs  in  the 
exact  expression  for  the  quantum  mechanical  transition  amplitude,  can  be 
neglected  because  each  separate  contribution  is  almost  always  dominant 
in  a temperature  region  in  which  the  other  contribution  is  negligible. 


III.  The  Overall  Rec-  ination  Coefficients  For  Low-Density  Plasmas 

In  the  zero-density  limit  the  overall  recombination  rate  coefficient 
simply  equals  the  total  rate  coefficient  for  all  spontaneous  radiative 
captures,  by  both  single-and  by  two-  electron  transitions 


i 


I 

i 
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a(Z,  i -♦  Z-l)  = a(Z,  i -»  Z-l;  k,  at) . (8) 

k,al 

At  any  given  finite  density  there  will  be  some  limiting  n-value  above 
which  subsequent  collisional  excitation  of  the  final  state  k,  at  becomes 
as  probable  as  radiative  decay  to  lower  levels.  It  can  be  argued  that 
levels  above  this  limit  should  not  be  included  in  equation  (8)s  but 
termination  of  the  summation  over  n at  the  limiting  value  obtained  from 
the  formula  given  by  Griem  (I9&O  probably  results  in  an  overestimate  of 
the  finite  density  reduction  of  the  recombination  rate.  In  the  present 
investigation  the  limiting  n-value  has  been  determined  by  the  less 
restrictive  condition  (Inglis  and  Teller  1939)  that  the  quasi-static 
linear-Stark  shifts  produced  by  all  surrounding  charged-particles  must 
not  exceed  the  separation  between  levels  with  adjacent  n- values. 

When  the  cut-off  procedure  affects  the  outer  electron  at-  levels 
which  play  the  most  important  role  in  the  dielectronic  recombination 
process,  a more  fundamental  treatment  of  the  density-dependence  is 
required.  Burgess  and  Summers  (1969)  solved  a set  of  statistical 
equilibrium  equations  which  describes  the  effects  of  collisions  on  the 
final  state  populations  and  calculated  a density-dependent  overall 
recombination  coefficient  by  a procedure  similar  to  that  employed  by 
Bates  et  al.  (1962) . However,  this  treatment  of  density-dependent 
dielectronic  recombination  is  incomplete  because  no  account  is  taken  of 
the  influence  of  the  plasma  on  the  doubly-excited  level  populations 
(Jacobs,  Davis  and  Kepple  1976) . 
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The  dielectronic  recombination  rate  coefficients  for  the  various 
ions  of  silicon  have  been  calculated  with  account  taken  of  all  stabilizing 
radiative  transitions  and  all  autoionization  processes  which  involve  a 
single-electron  electric-dipole  transition  of  the  recombining  ion  core. 

A distorted  wave  calculation  (Davis  et  al.  1977)  was  employed  to  obtain 
the  threshold  values  of  the  collisional  excitation  cross  sections, 
because  the  widely-used  Bethe  approximation  is  known  to  greatly  over- 
estimate the  most  important  partial  wave  contributions.  In  Table  I the 
stabilizing  radiative  transitions  j -»  k are  given  in  order  of  decreasing 
importance  at  the  temperature  where  the  ion  in  question  is  predicted  by 
the  present  calculations  to  have  its  maximum  equilibrium  abundance.  The 
asterisk  has  been  used  to  identify  the  transitions  which  involve  doubly- 
excited  states  subject  to  autoionization  into  an  excited  state  of  the 
recombining  ion. 

In  Table  II  the  dielectronic  contributions  to  the  overall  recom- 
bination rate  coefficients  are  presented  as  functions  of  log^Q  Te(°K) . 

Also  given  are  the  values  which  were  obtained  by  evaluating  the  Burgess 
formula  with  our  nonhydrogenic  energy  eigenvalues  and  oscillator  strengths 
for  the  transitions  in  Table  I.  The  effects  of  the  surrounding  charged 
particles  in  a plasma  with  an  electron  density  of  10*^  cm  ^ have  been 
estimated  by  terminating  the  summation  over  n in  equation  (8)  at  the 

Inglis-Teller  limit.  The  limiting  n-value  obtained  from  the  Inglis- 

3/5 

Teller  criterion  increases  with  Z as  Z , whereas  for  a given  stabilizing 


radiative  transtion  j k the  n- values  which  give  the  dominant  contri- 
bution to  the  total  dielectronic  recombination  rate  are  descreasing  with 
increasing  Z.  Consequently,  the  recombination  rates  obtained  for  the  most 


r 
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highly-charged  ions  are  to  a good  approximation  independent  of  density 
over  an  extensive  density  range.  Three-body  recombination  is  negligible 
at  the  densities  of  interest. 

In  the  temperature  regions  of  maximum  equilibrium  abundance,  the 
dielectronic  recombination  rate  coefficients  obtained  from  the  present 
calculation  are  in  several  cases  substantially  reduced  in  comparison 
with  the  values  obtained  from  the  Burgess  formula.  For  Si  V - Si  X 
this  reduction  is  attributable  to  the  effect  of  autoionization  into  the 
3p  excited  state  of  the  recombining  ion,  which  is  found  to  be  the 
dominant  autoionization  process  for  the  3d,  n-t  doubly-excited  states. 

The  effect  of  autoionization  accompanied  by  excitation  is  illustrated 
in  Figure  2 by  showing  as  functions  of  log10  Te(°K)  the  contributions  to 
the  Si  VI  overall  recombination  coefficient  arising  from  the  direct  radi- 
ative transition  and  from  the  three  most  important  dielectronic  transitions. 
When  no  account  is  taken  of  autoionization  into  the  3p  oxci ted  state,  the 
dielectronic  contribution  arising  from  the  3d  ->  2p  radiative  transition 
is  found  to  be  dominant  in  the  maximum  equilibrium  abundance  temperature 
region.  This  is  illustrated  by  the  values  obtained  using  the  Burgess 
formula,  which  are  in  this  case  about  two  orders  of  magnitude  larger  than 
the  present  results.  With  the  inclusion  of  the  additional  3d  -*  3p  auto- 
ionization rate,  the  3s  -*  2p  and  2p  -*  2s  dielectronic  contributions  are 
relatively  more  important. 


IV.  Corona  Ionization  Equilibrium 

The  corona  equilibrium  distributions  of  the  ions  of  a given  element 


Li 


are  determined  by  the  conditions 
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N(Z-l,g)  S(Z-l,g  -+  Z)  = N (Z,g)  a(Z,g-*Z-l),  UZiZ  (9) 

max 

together  with  the  normalization  requirement 
^max 

NT  = S NCZ,g>  (10) 

Z = 0 

N(Z,g)  is  the  density  of  ions  with  residual  charge  Z in  the  ground  state, 

while  Nt  is  the  total  density.  S(Z-1,  g -»  Z)  denotes  the  total  rate 

(Z-l) 

coefficient  for  single  ionization  from  the  ground  state  of  X as  a 

result  of  thermal  electron  impacts.  Q!(Z,  g -4  Z-l)  is  the  overall  recom- 

(Z) 

bination  rate  coefficient  from  the  ground  state  of  X ' . When  three-body 
recombination  and  the  density  dependence  of  the  spontaneous  radiative 
capture  rate  discussed  in  the  preceding  section  are  neglected,  the 
relative  abundances  N (Z,g)/NT  are  independent  of  density  and  are 
functions  only  of  the  local  electron  temperature. 

The  ionization  rate  coefficients  employed  in  the  present  calculations 
were  obtained  by  evaluating  the  semi-empirical  expression  given  by 
Lotz  (I967).  Contributions  from  inner-shell  ionization  have  been  included. 
The  additional  contributions  from  autoionization  following  inner-shell 
excitation  have  been  estimated  by  using  a procedure  similar  to  the  one 
employed  by  Jordan  (I969).  Experimentally  determined  (Datla,  Nugent,  and 
Griem  1976)  ionization  rate  coefficients  for  certain  multiply-charged 
ions  are  as  much  as  a factor  of  2 less  than  the  predictions  of  the  semi- 
empirical  formula.  This  uncertainty  is,  however,  small  compared  to  the 
reductions  which  have  been  obtained  in  the  overall  recombination  rates 
by  the  inclusion  of  autoionization  accompanied  by  excitation. 


MM . 
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The  relative  steady-state  abundances  of  the  various  ions  of  silicon 
are  presented  in  Table  III  and  shown  in  Figure  3 as  functions  of 
log^  (°K) . The  maximum  equilibrium  abundance  temperatures  are  given 
in  Table  IV  together  with  those  predicted  by  the  calculations  of  Jordan 
(I969)  and  Summers  (197*0.  The  differences  between  the  results  of  the 
two  previous  calculations  are  most  probably  attributable  to  the  use  of 
different  ionization  rate  coefficients.  Very  close  agreement  with  the 
results  of  Jordan  (1969)  is  obtained  in  the  present  calculations  by  the 
omission  of  the  rates  for  autoionization  into  an  excited  state.  Conse- 
quently the  shifts  in  the  abundance  curves  toward  lower  temperatures, 
which  are  illustrated  for  Si  V - Si  VII  in  Figure  4,  are  directly 
attributable  to  the  effects  of  autoionization  accompanied  by  excitation. 
The  reduced  maximum  abundance  temperatures  which  have  been  obtained  for 
certain  ions  are  still  significantly  higher  than  those  predicted  (House 
1964)  when  no  account  is  taken  of  the  effects  of  dielectronic  recombi- 
nation. It  should  be  emphasized  that  a relatively  small  shift  of  a 
sharply  peaked  abundance  curve  can  significantly  alter  the  predicted 
emission  line  intensity. 

V.  Radiative  Energy  Loss  Rates 

In  this  section  the  results  are  presented  for  the  rates  of  radiative 
energy  loss  arising  from  radiation  processes  which  involve  collisions 
between  silicon  ions  and  thermal  electrons.  The  ions  are  assumed  to  be 
in  corona  ionization  equilibrium,  and  no  account  is  taken  of  the  effects 
of  self-absorption. 
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The  rate  at  which  energy  is  radiated  per  unit  volume  as  a result  of 
the  electron  impact  excitation  of  the  resonance  lines  is  given  by 

Z -1 
max 

€L  = Ne  ^ Z AE(Z,g  -»  j)  C(Z,g  -»  j)  B(Z,j  -»  g)  (11) 

z = l Nt  J 

The  excitation  energies  are  denoted  by  AE(Z,g  -»  j),  and  C(Z,g  -*  j)  are 

* 

the  collisional  excitation  rate  coefficients  obtained  from  the  distorted 
wave  calculation  (Davis  et  al.  1977)-  The  branching  ratios  B(z,j  -*  g) 
for  the  resonance  line  emissions  are  very  close  to  unity. 

The  total  power  radiated  per  unit  volume  during  dielectronic  recom- 
bination is  given  by 

Z -1 
max 

ed  = \ Nt  Z Z AE(z,g  -♦  j)  ad(z’g;  j -»  g),  (12) 

Z = 1 nt  j 

where  (Z,  g;  j -♦  g)  denotes  the  contribution  to  the  overall  dielec- 
tronic recombination  rate  coefficient  arising  from  the  stabilizing  radi- 

I 

ative  transition  j g.  No  account  has  been  taken  of  the  radiation  which 
is  emitted  when  the  nA-states  of  the  outer  electron  cascade  to  the  ground 

| | 

states.  In  addition,  we  have  neglected  the  displacement  of  the  emitted 

; 

photon  frequency  from  the  center  of  the  resonance  line  which  results  from 

! 

the  presence  of  the  outer  electron.  The  single-electron  electric-dipole 
transitions  in  Table  I have  been  taken  into  account  in  the  evaluation  of 
both  the  resonance  line  and  the  dielectronic  recombination  contributions. 

The  results  are  presented  in  Figure  5 By  plotting  as  functions  of 


loglrt  T (°K)  the  rate  coefficients  /N  Nm  and  £,/N  N_,  which  have  the 

iUe  LeTdeT 

dimensions  of  power  X volume.  A^so  shown  are  the  corresponding  rate 
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coefficients  for  direct  recombination  radiation  and  for  bremsstrahlung 
which  were  evaluated  as  described  by  Griem  (196^).  The  relative  maximum 
produced  by  the  2p  -*  2s  transitions  in  Si  VI  - Si  XII  is  found  to  occur 
very  close  to  106  °K,  which  is  lower  than  predicted  by  the  calculations 
of  Cox  and  Tucker  (1969).  This  reduction  is  directly  attributable  to 
the  effects  of  the  reduced  recombination  rates  on  the  ionization 
equilibrium  distributions. 

In  Table  V the  individual  contributions  to  the  rate  coefficients 
^l/^e  NT  an<*  ^d^Ne  ^T  ^rom  most  important  radiative  transitions  j -»  g 
are  given  at  the  maximum  line  emission  temperatures,  which  are  in  several 
cases  significantly  higher  than  the  maximum  abundance  temperatures  in 
Table  IV.  The  intensity  of  dielectronic  recombination  radiation,  relative 
to  the  associated  resonance  line  intensity,  is  found  to  be  greatest  in 
the  cases  of  the  2p  -»  Is  transitions  of  the  helium-like  and  hydrogen-like 
ions  Si  XIII  and  Si  XIV.  It  has  been  pointed  out  by  Gabriel  (I972)  that 
for  higher  members  of  the  helium  isoelectronic  sequence  the  dielectronic 
recombination  process  occurs  primarily  through  the  n = 2 levels  of  the 
outer  electron,  which  give  rise  to  satellite  features  that  are  spectros- 
copically resolvable  from  the  resonance  line. 

Although  the  radiation  emitted  during  dielectronic  recombination  can 
be  considerably  more  important  than  direct  recombination  radiation  and 
bremsstrahlung,  the  electron  impact  excitation  of  resonance  line  radiation 
is  clearly  the  dominant  radiative  cooling  mechanism  at  temperatures  where 
ions  with  bound  electrons  are  abundant.  In  a rapidly  ionizing  plasma  the 
ionization  of  multiply-charged  ions  may  represent  an  important  energy 
loss  mechanism.  Resonance  line  emission  is  expected  to  be  the  dominant 
energy  loss  mechanism  for  multiply-charged  ions  in  ionization  equilibrium. 
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VI.  Conclusions 

With  account  taken  of  the  effects  of  autoionization  accompanied  by 
excitation,  the  present  calculation  for  the  dielectronic  recombination 
rates  and  the  ionization  equilibrium  distributions  may  represent  the 
best  improvement  that  is  possible  within  the  single-electron  electric- 
dipole  and  simple  corona  model  approximations.  The  use  of  a more 
sophisticated  calculation  for  the  ionization  rates  would  be  desirable. 
However,  the  most  significant  additional  improvements  are  expected  to 
result  from  taking  into  account  more  general  autoionization  and  radiation 
procesess  in  a treatment  of  the  effects  of  collisions  on  both  the  singly- 
and  the  doubly-  excited  level  populations.  This  more  sophisticated 
treatment  of  dielectronic  recombination  will  be  the  subject  of  a future 
investigation. 
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Table  I 


Stabilizing  Radiative  Transitions  in  Order  of  Decreasing 
Importance  at  the  Temperature  of  Maximum  Equilibrium  Abundance 


[Recombining 
Ion 


Single  - Electron  Transitions 


Si 

II 

3p 

-4 

3s, 

4d 

-4 

3p*, 

bs 

-4 

3p 

Si 

III 

3p 

-4 

3s, 

4p 

-4 

3s* 

SI 

IV 

3p 

-4 

3s, 

4p 

-4 

3s* 

Si 

V 

3s 

-4 

2p, 

3d 

-4 

2p* 

Si 

VI 

3s 

-4 

2P, 

2p 

-4 

2s, 

3d 

-4 

2p* 

Si 

VII 

2p 

-4 

2s, 

3d 

-4 

2p*» 

3s 

-4 

2p 

Si 

VIII 

2p 

-4 

2s, 

3d 

-4 

2p* , 

3s 

-4 

2p 

Si 

IX 

2p 

-4 

2s, 

3d 

-4 

2p* , 

3s 

-4 

2p 

Si 

X 

2p 

-4 

2s, 

3d 

-4 

2p*, 

3s 

-4 

2p 

Si 

XI 

2p 

-4 

2s, 

3p 

-4 

2s* 

Si 

XII 

2p 

-4 

2s, 

3p 

-4 

2s* 

Si 

XIII 

2p 

-4 

Is* 

3d  -»  3p 


t Affected  by  Autoionization  Accompanied  by  Excitation 
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Table  II 

Dielectronic  Recombination  Rate  Coefficients  (Cm3  sec-1) 


Recombining  Ions 


Si  V 


Si  II 

Si  III 

Si  IV 

o.i6(-ii)* 

o.44(-io)** 

0.49(-12) 

o.66(-li) 

0.54(-U) 

0.56(-12) 

o.54(-ii) 

o.77(-io) 

0.48(-ll) 

0.46(-l0) 

O.ll(-lO) 

0.12(-10) 

0.34(-l0) 

o.69(-10) 

0.66(-lo) 

0.16C-9) 

o.i3(-io) 

0.99C-10) 

0.27(-10) 

0.l4(-9) 

0.20(-10) 

0.10(-9) 

o.70(-io) 

o.26(-9) 

0.20(-10) 

0.84(-10) 

0.10(-9) 

0.32(-9) 

0.81(-10) 

0.22C-9) 

o.i6(-io) 

o.59(-io) 

0.99(-l0) 

0.28C-9) 

0.72(-10) 

0.21C-9) 

O.ll(-lO) 

0.37(-l0) 

o.76(-io) 

0.20C-9) 

0.52(-10) 

0.l6(-9) 

0.65(-ll) 

0.22(-10) 

o.50(-io) 

0.l3(-9) 

0.33C-10) 

0.1K-9) 

0.36(-ii) 

0.12(-10) 

0.30C-10) 

o.76(-io) 

o.i9(-io) 

o.63(-io) 

0.20(-ll) 

0.64(-ll) 

o.17(-io) 

0.42(-10) 

o.io(-io) 

0.35("10) 

O.lO(-ll) 

o.33(-n) 

o.89(-ii) 

0.22(-10) 

0.55(-ll) 

0.19(-10) 

o.53(-l2) 

0.17(-H) 

0.47(-ll) 

0.12(-10) 

0.29(-ll) 

0.98(-ll) 

0.27(-12) 

0.87(-12) 

0.24(-ll) 

o.6o(-ii) 

0.15(-U) 

0.50C-11) 

0.14(-12) 

0.44C-12) 

0.12(-11) 

0.30C-11) 

o.75(-12) 

0.26(-ll) 

.63(-37) |0.15(-17) 
0.42(-33) 

0.17C-15) 
0.17 (-24) 

0.26(-l4) 
0.36C-19) 

0.12C-13) 

0.65C-16) 

0.30C-13) 

o.6o(-i4) 

0. 10(-12) 
0.14(-12) 

o.34(-12) 
0.19C-H) 

0.84C-12) 

o.io(-io) 

0.17C-H) 

0.28(-10) 

o.27(-ii) 
0.45C-10) 

0.29C-11) 
0.49C-10) 

o.24(-ii) 

o.4i(-io) 

0.17C-H) 

0.28C-10) 

O.lO(-ll) 
0.18(-10) 


Si  VII 


0.18C-19) 

0.33(-3D 

0.44C-16) 

0.10(-22) 

0.47C-14) 

0.18C-17) 

0.76C-13) 
0.29 (-14) 

0.56(-l2) 

0.23C-12) 

o.34(-ii) 

0.29(-ll) 

o.ii(-io) 

o.ii(-io) 

0.20(-10) 

0.23C-10) 

0.23C-10) 

o.4o(-io) 

0.21 (-10) 
o.6o(-io) 

0.17(-10) 

0.68(-10) 

0.12(-10) 

0.59(-10) 

o.76(-11) 

0.43(-lO) 

0.45(-ll) 

0.27(-10) 


Table  II  (continued 


°g,«T  (°K) 


5.0 


Recombining  Ions 

Si  VIII 

Si  IX 

Si  X 

Si  XI 

Si  XII 

Si  XIII 

0.13C-10) 

0.12C-10) 

o.27(-io) 

o.23(-io) 

o.37(-io) 

0.34(-l0) 

0.43 (-10) 
0.33C-10) 

0.36(-l0) 

0.47C-10) 

0.27C-10) 

0.28C-10) 

o.46(-lo) 

0.47C-10) 

o.65(-lo) 

0.68(-10) 

0.78(-10) 

0.78(-i0) 

o.44(-io) 

0.57C-10) 

0.36C-10) 

0.38C-10) 

o.55(-io) 

0.58C-10) 

0.77C-10) 

0.82(-10) 

0.10C-9) 

0.10(-9) 

o.4o(-io) 

0.51C-10) 

O.lk(-lO) 

o.43(-io) 

o.49(-io) 

o.54(-io) 

0.68(-l0) 

0.72(-10) 

0.94(-l0) 

0.98(-10) 

0.29(-10) 

0.37C-10) 

0.79 (-28) 

0.45C-29) 

0.28(-l0) 

0.53C-10) 

0.38C-10) 

0.53C-10) 

o.50(-io) 

0.58(-10) 

0.7H-10) 

0.77C-10) 

0.19(-10) 

0.25(-l0) 

0. 11(-21) 
0.82 (-22) 

0.21(-10) 

0.63C-10) 

0.27(-10) 

0.57(-lO) 

0.3'+(-l0) 

0.48(-10) 

0.46(-lo) 

o.6o(-io) 

O.ll(-lO) 

0.21(-10) 

0.7K-18) 

o.24(-l7) 

0.15C-10) 

o.59(-io) 

0.19C-10) 

0.52(-10) 

0.22(-l0) 

o.4o(-io) 

0.28(-10) 

o.49(-io) 

o.69(-io) 

0.20(-10) 

0.2K-15) 

0.13C-14) 

o.ioc-io) 

o.46(-io) 

0.13C-10) 

0.42(-10) 

0.l4(-lo) 

0.31(-10) 

o.i6(-io) 

0.37C-10) 

0.42(-ll) 

0.18(-10) 

0.85(-l4) 

0.52C-13) 

o.6i(-H) 

0.30(-10) 

0.78C-11) 

0.28(-10) 

0.84(-ll) 

0.21(-10) 

0.86(-ll) 

0.25(-10) 

0.24(-ll) 

0.13(-10) 

o.82(-l3) 

o.43(-l2) 

o.35(-H) 

0.18(-10) 

0.4-5C-11) 

0.17(-10) 

0.48(-ll) 

o.i3(-lo) 

o.46(-H) 

o.i6(-io) 

0.14(-11) 

0.87C-11) 

0.28(-12) 

0.12(-11) 

o.i9(-ii) 

0.10(-10) 

0.25 (-11) 

0.98(-ll) 

0.26(-ll) 

0.72(-ll) 

0.24(-ll) 

0.92(-ll) 

o.76(-12) 

0.52(-ll) 

0.49(-12) 

0.19(-H) 

O.lO(-ll) 

o.55(-n) 

0.13(-U) 

0.53(-ll) 

o.l4(-il) 

0.39(-ll) 

0.12(-11) 

0.50(-ll) 

0.40(-12) 

o.29(-ii) 

0.52(-12) 

0.19C-H) 

0.52C-12) 

0.29C-11) 

0.69C-12) 

0.28(-ll) 

0.71(-12) 

0.21(-11) 

0.63(-12) 

0.27(-ll) 

0.21(-12) 

0.l6(-ll) 

0.44(-12) 

0.15(-H) 

0.27C-12) 

o.i5(-ii) 

0.35C-12) 

o.i4(-li) 

0.37C-12) 

O.ll(-ll) 

0.32(-12) 

0.l4(-ll) 

0. 11(-12) 
0.82(-12) 

o.30(-12) 

0.99C-12) 

Numbers  in  parentheses  are  the  powers  of  10. 
Calculated  using  the  formula  of  Burgess  (I965) 


woe********  ■ 
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Table  III 

-log^0  N(Z)/Nt  for  Silicon  Ions 


^loV0® 

Si  II 

Si  III 

Si  IV 

Si  V 

4.0 

0.06 

4.34 

4.1 

0.02 

2.76 

4.2 

0.01 

1.73 

4.3 

0.04 

1.04 

4.4 

0.13 

0.58 

5.00 

4.5 

0.30 

0.29 

3.43 

4.6 

0.56 

0.14 

2.25 

4.94 

4.7 

0.87 

0.09 

1.34 

2.74 

4.8 

1.28 

0.17 

0.72 

1.05 

4.9 

2.20 

0.77 

0.73 

0.19 

5.0 

3.54 

1.84 

1.27 

0.03 

5.1 

4.82 

2.86 

1.85 

0.01 

5.2 

3.73 

2.33 

0.01 

los10T.<=K> 

Si  V 

Si  VI 

Si  VII 

Si  VIII 

Si  IX 

5.0 

0.03 

5.1 

0.01 

5.2 

0.01 

5.3 

o.o4 

4.30 

5.4 

0.23 

0.37 

2.58 

5.5 

0.66 

0.12 

1.49 

3.93 

5.6 

1.17 

0.11 

0.77 

2.34 

5.00 

5.7 

1.76 

0.30 

0.37 

1.21 

2.98 

5.8 

2.53 

0.72 

0.30 

0.54 

1.57 

5.9 

3.50 

1.41 

0.57 

0.31 

0.73 

6.0 

4.70 

2.34 

1.14 

0.45 

0.37 

6.1 

3.51 

1.99 

0.95 

0.45 

6.2 

5.00 

3.24 

1.88 

1.02 
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Table  IV 

Maximum  Equilibrium  Abundance  Temperatures (in  °K)for  Silicon  Ions 


i 

Ion 

Present  Results 

Jordan  (I969) 

Summers  (1974) 

Si 

II 

1.6(4) 

1.6(4) 

1.6(4) 

Si 

III 

5.0(4) 

5.0(4) 

6.3(4) 

Si 

IV 

6.3(4) 

8.0(4) 

8.0(4) 

Si 

V 

1.6(5) 

2.0(5) 

2.0(5) 

Si 

VI 

4.0(5) 

5.0(5) 

5.0(5) 

Si 

VII 

6.3(5) 

7.0(5) 

1.0(6) 

Si 

VIII 

8.0(5) 

1.0(6) 

1.3(6) 

Si 

IX 

1.0(6) 

1.3(6) 

1 .6(6) 

Si 

X 

1.3(6) 

1 .6(6) 

2.0(6) 

Si 

XI 

1.6(6) 

2.0(6) 

2.5(6) 

Si 

XII 

2.0(6) 

2.0(6) 

2.5(6) 

Si 

XIII 

4.0(6) 

4.0(6) 

4.0(6) 

Si 

XIV 

1.3(6) 

1.3(7) 

1.3(7) 

Table  V 

Individual  Contributions  to  the  Energy  Loss  Coefficients  for 
Dielectronic  Recombination  and  Line  Radiation,  evaluated  at 
the  Maximum  Line  Emission  Temperature 


Ion 

Radiative  Transition 

J -k# 

Si  II 

3p  -♦  3s 

4d  -♦  3p 

Si  III 

3p  -♦  3s 

4p  -»  3s 

Si  IV 

3p  -*  3s 

4p  -+  3s 

i 

Si  V 

3s  -»  2p 

3d  -»  2p 

t Si  VI 

2p  -»  2s 

3s  -*  2p 

3d  -»  2p 

Maximum 
Line  Emission 
Temperature 
°K 


6.3(4) 

6.3(4) 


8.0(4) 

8.0(4) 


Wt  €L/Me\ 
(Watts  cm3)  (Watts  cm3) 


0.90(-29) 

0.58C-30) 


0. IK-27) 

0.80(-30) 


0.32(-28) 

0.20(-30) 


0.49(-29) 

0.34(-30) 


0.36(-29) 
0. ll(-28) 
0.49(-29) 


0.24 (-25) 
0.12(-25) 


0.12(-24) 

0.12(-24) 


0.22(-25) 

0.96(-26) 


0.1o(-27) 

0.57C-27) 


0.33(-27) 

0.50(-27) 

0.36(-26) 


-24- 

Table  V (continued 


Si  VII 

2p  -*  2s 

6.3(5) 

o.75(-28) 

0.90(-2 6) 

3s  -»  2p 

6.3(5) 

0.13(-28) 

o.46(-27) 

3d  -»  2p 

6.3(5) 

0.23 (-28) 

0.55(-26) 

Si  VIII 

2p  -»  2s 

8.0(5) 

0.60(~28) 

0.l3(-25) 

3s  -»  2p 

1.0(6) 

0.12(-28) 

0.45(-27) 

3d  -»  2p 

1.0(6) 

0.42(-28) 

0.63 (-26) 

Si  IX 

2p  -*  2s 

1.0(6) 

0.54(-28) 

o.l3(-25) 

3s  -»  2p 

1.3(6) 

0.97(-29) 

0.32(-27) 

3d  -►  2p 

1.3(6) 

0.64(-28) 

0.57(-26) 

Si  X 

2p  -»  2s 

1.3(6) 

0.47(-28) 

0.l4(-25) 

3s  -»  2p 

1.6(6) 

0.40(-29) 

o.  IK-27) 

3d  -»  2p 

1.6(6) 

0.42(-28) 

0.27(-26) 

Si  XI 

2p  -♦  2s 

1.6(6) 

0.45 (-28) 

0.12(-25) 

3p  -»  2s 

l .6(6) 

0.85(-29) 

0.20(-26) 

Si  XII 

2p  -*  2s 

2.0(6) 

0.35(-29) 

o.38(-26) 

3p  -»  2s 

2.0(6) 

0.13(-28) 

0.14 (-26) 

Si  XIII 

2p  -»  Is 

1.0(6) 

0.66(-28) 

0.22(-26) 

Si  XIV 

2p  -»  Is 

1.6(6) 

0.i6(-27) 

0.13(-26) 

* Including  all  LS-  Components. 

•Numbers  in  parentheses  are  the  powers  of  10. 
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Figure  Captions 

Figure  1.  Radiationless  capture,  autoionization,  and  radiative  decay 
processes  determining  the  3d,  nt  doubly-excited  level 
populations  in  boron-like  through  neon-like  ions. 

Figure  2.  Direct  radiative  and  dielectronic  contributions  to  the 
Si  VI  recombination  rate  coefficient. 

present  calculation 

___________  Burgess  formula  (I965) . 

The  corona  ionization  equilibrium  of  silicon  ions. 

The  corona  ionization  equilibrium  distributions  for 
Si  V - Si  VII. 

present  calculation 

-----------  Jordan  (I969) . 

Figure  J.  The  radiative  energy  loss  rate  coefficients  for  silicon  ions 
in  isothermal  steady-state  plasmas. 
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